Dynamical Scaling: the Two-Dimensional XY Model Following a Quench 
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To sensitively test scaling in the 2D XY model quenched from high-temperatures into the or- 
dered phase, we study the difference between measured correlations and the (scaling) results of a 
Gaussian-closure approximation. We also directly compare various length-scales. All of our results 
are consistent with dynamical scaling and an asymptotic growth law L ~ (t/ln[t/to])^ , though 
with a time-scale io that depends on the length-scale in question. We then reconstruct correlations 
from the minimal-energy configuration consistent with the vortex positions, and find them signifi- 
cantly different from the "natural" correlations — though both scale with L. This indicates that 
both topological (vortex) and non-topological ( "spin- wave" ) contributions to correlations are rele- 
vant arbitrarily late after the quench. We also present a consistent definition of dynamical scaling 
applicable more generally, and emphasize how to generalize our approach to other quenched sys- 
tems where dynamical scaling is in question. Our approach directly applies to planar liquid-crystal 
systems. 
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I. INTRODUCTION 

The study of non-equilibrium dynamics in systems 
with continuous symmetries has burgeoned |l|. Liquid- 
crystalline systems [§|-^, evolving after being quenched 
into an ordered phase, provide picturesque examples of 
topological defects and their interactions. Evolving sys- 
tems of topological defects arc also found in applications 
from cosmology ||] to quantum Hall ferromagnets ||l^,|ll| . 

A relatively simple system with a continuous symmetry 
is the two-dimensional XY ferromagnet with no disorder, 
which supports singular vortices that carry topological 
charge and have logarithmic interactions. The equilib- 
rium properties have spawned a rich and fertile litera- 
ture punctuated by the work of Kosterlitz and Thou- 
less [ |l2[ . More recently, the non-equilibrium behavior 
of the 2D XY model following a quench to below the 
Kosterlitz- Thouless critical temperature, Tkt, has been 
studied theoretically p"3|-p0{ and also experimentally [plH] 
with specially prepared liquid-crystal systems. Related 
2D liquid-crystal systems have also been studied theoret- 
ically pl]-|23|| and experimentally [^,HJ|]. 

Following a quench at i = from a disordered phase 
into an ordered phase, a crucial issue is whether there is 
dynamical scaling |24] at late times t, where 



Cir, t) = {4>{x, t) ■ 0(x -I- r, t)) = f{r/L). 



(1) 



Here, (p is the XY order parameter, /(x) is a time- 
independent scaling-function for the two-point correla- 
tions, and L{t) is a growing length-scale that captures 



all of the correlation dynamics. The explicit or implicit 
assumption of dynamical scaling underpins most theo- 
retical descriptions of phase-ordering structure p|,p5[-p7[ . 
Unfortunately, apart from a limited number of solvable 
systems, there exist no theoretical approaches to a pri- 
ori determine dynamical scaling. Indeed, the presence 
or absence of dynamical-scaling remains an unresolved 
issue in the 2D XY model |1^J2^. This is surprising, 
since simple systems that break scaling are seen as ex- 
ceptions p8| . For example, the weak scaling violations in 
the conserved spherical model identified by Coniglio and 
Zannetti p9[ are due to non-commuting spherical and 
asymptotic-time limits 1 30 related to similar phenomena 
in equilibrium critical dynamics [ pl| . 

Stronger scaling violations are found in one- and two- 
dimensional systems with non-singular topological tex- 
tures [ pT|j32| . These systems segregate into domains of 
similarly charged textures, similar to the morphologies 
seen in reaction-diffusion A + B ^> % systems p3[ . The 
domain-size and the texture separation provide distinct 
growing length-scales. Within this context, the difficulty 
in resolving scaling in the 2D XY model can be under- 
stood. Viewed as a plasma of overdamped charged vor- 
tices with logarithmic interactions [p4| , quenched from 
high-temperatures, the 2D XY model sits exactly at the 
marginal dimension [d = 2) below which segregated mor- 
phologies with strong scaling violations are expected, and 
above which a mixed morphology with only one length- 
scale, the particle separation, is seen [nOl. Such particle 
systems are expected to scale, with no domain structure. 
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at the marginal dimension |L0| , however the asymptotic 
regime could be quite late. 

With dissipative dynamics and the assumption of 
dynamical-scaling the predicted as ymp totic growth-law 
of the characteristic length-scale is 



L{t)^A{t/Ht/to]Y'\ 



(2) 



where A and io are the non-universal amplitude and time- 
scale, respectively. This growth-law characterizes the 
correlations with a length Li/2(i), where C[Li/2,i) = 
1/2, as well as the vortex separation with a length-scale 
Ly{t), where the vortex density pdef — l/^u- These 
lengths will only differ by prefactors and by subdominant 
contributions at late times. [Eqn. also describes the an- 
nihilation time of an isolated vortex-antivortex pair with 
an initial separation L ||l^.] The logarithmic factor is 
crucial, and stems from the logarithmic vortex mobility. 
The same growth-law is expected in liquid-crystal films 
with vortices |^ . 

The analytical evidence for scaling violations is mostly 
suggestive: explicit violations in four-point correlations 
pM and multiple energy-scales seen in energy-scaling cal- 
culations [g3- These would indicate multiple lengths 
which differ at most by logarithmic factors, consis- 
tent with the marginal dimensionality within a reaction- 
diffusion context [nol. Indeed, approximation schemes for 
correlation functions in the 2D XY model typically find 
scaling but with no logarithmic factors (see, e.g., P6| , p5| , 
see also |Q). Additionally, the 2D XY model quenched 
between two temperatures below Tkt, and coarse-grained 
to a fixed scale to eliminate bound vortex pairs, is solv- 
able |3^ and dynamically scales without any logarithmic 
factor, L{t) - t^/"^. 

Previous numerical evidence for scaling violations is 
stronger. Cell-dynamical simulations of XY models 
quenched to T = by Blundell and Bray |l^ found that 
two-point correlations did not scale well with respect to 
the defect separation Ly, though they scale with respect 
to the correlation length L1/2 (see also |13pG]). Mon- 
dello and Goldenfeld ||l^ also found indications of multi- 
ple length-scales. Simulations of nematic films by Zapo- 
tocky et al. [g2| found a variety of effective growth expo- 
nents, though again the correlation function appeared to 
scale (see also pl| , p3| ). Other simulations on the 2D XY 
model at finite temperatures have recovered the expected 
growth law |l7||l9| ], and have found dynamical scaling 
||19| . Simulations of quenches to T = in hard-spin sys- 
tems found dynamical scaling of correlations even though 
the dynamics froze at late times |14[| ! 

Experiments on liquid-crystal systems, following the 
pioneering work by Shiwaku et al. [^ , have recovered the 
^1/2 growth of defect separation after a quench, though 
with insufficient resolution to determine logarithmic fac- 
tors 0Tj,Hl and with some difficulties in achieving an 
unbiased (symmetric) quench nM- When measured, the 
structure [^pj and other two-point correlations are con- 
sistent with dynamical scaling 0] . 



In this paper we want to clarify the existence or ab- 
sence of dynamical scaling in the 2D XY model. A suc- 
cessful strategy can then be applied more generally to 
systems that seem to violate scaling, in particular to sys- 
tems with more complicated collections of defects p8] . 

We first discuss the appropriate definition of dynami- 
cal scaling, within the context of systems relaxing after 
a quench. We then derive approximate forms for vari- 
ous correlation functions via Gaussian closure techniques, 
which impose scaling. While we do not expect them to 
exactly match the measured correlations, they are used 
to normalize the measured values in order to enhance our 
sensitivity to scaling or its absence. In combination with 
the growth-law, we have a "null hypothesis" which would 
be broken by scaling violations. We present our simula- 
tion data and find no evidence for scaling violations. We 
then explicitly reconstruct a two-point gradient correla- 
tion function, within the periodic system using only the 
vortex positions and charges, and find it significantly dif- 
ferent from the unreconstructed scaling form. However 
both correlations scale with respect to the defect den- 
sity. This indicates that both topological (vortex) and 
non-topological ( "spin- wave" ) contributions to the order 
parameter are asymptotically relevant, with characteris- 
tic lengths that remain asymptotically proportional. 



II. DYNAMICAL SCALING 



In phase-ordering, dynamical scaling colloquially 
means that there is a single characteristic length scale 
growing in time. This leads to a rough-and-ready symp- 
tom of dynamical scaling violations: multiple length- 
scales with distinct growth-laws, see for example p2| , p8[ . 
While useful as a guide, this approach has limitations. 
One must first identify each asymptotic growth law, i.e. 
the effective exponent after it is constant in time and 
before finite-size effects of the sample become impor- 
tant. Practically, at most one or two decades in time 
are available in simulations if a 5% exponent variation is 
tolerated, and often less than a decade in experiments. 
When the scaling prediction for the growth-law is not a 
priori known, this approach on its own is dangerous. In- 
deed, sub-dominant corrections to the asymptotic growth 
law p^ can depend on the method used to extract the 
length-scale |Q . Even the observation of two asymptot- 
ically distinct length-scales does not demonstrate that 
they are dynamically interconnected. A silly example 
helps here: consider a sample made from gluing together 
a conserved binary-alloy system (asymptotic growth law 
t^'^), and a non-conserved order-disorder alloy system 
(growth law t^/^). Clearly two-growth laws could be ob- 
served in the hybrid, but they should not imply scaling 
violations. [Such dynamically independent sub-systems 
would lead to correlation functions that are sums of scal- 
ing functions.] The situation is more complicated when 
both lengths are observed within a homogeneous sample. 



such as the asymptotic behavior of monopoles and vortex 
lines in bulk nematics [|6|. Non-trivial inter-relationships 
of observed lengths can generally only be resolved with 
the help of simplified dynamical models, for example see 

A more precise definition of dynamical scaling is 
that two-point equal-time correlations have a time- 
independent scaling form, see Eqn. H^, which also implies 
scaling of the structure factor 



S{k, t) = {${k, t) ■ ${~k, t)) = L'^gikL), 



(3) 



where g{x) is a time-independent scaling function. This 
is directly measured in scattering experiments, can be 
well approximated analytically, and is easy to extract 
from simulations. For systems with singular topologi- 
cal defects, such as domain walls, hedgehogs, vortices, 
or vortex lines, a generalized Porods law 0| connects the 
density of defect core pdef to the asymptotics of the struc- 
ture via 



5(fc) ^ Pde/A:"('^+"Ui » 1, 



(4) 



where n characterizes the defect type [for the 2D XY 
model, n = d = 2]. This directly implies that the length 
derived from the defect density, Ly , is asymptotically pro- 
portional to the correlation length L1/2 when the corre- 
lations dynamically scale. 

This definition is still incomplete, since systems can 
satisfy Eqn. yet have distinct lengths intimately con- 
nected by the dynamics — e.g. in the ID XY model 
[ p2[ . Additionally, higher-point correlations can be con- 
structed in the 2D XY model which explicitly do not scale 
p6| , p2[ . Should these be viewed as violations of dynam- 
ical scaling? Fortunately a self-contained definition of 
dynamical scaling exists, introduced by Bray and Rutcn- 
berg [2^ . In order to calculate the rate of free-energy dis- 
sipation in a coarsening system, they additionally require 
the scaling of the time-derivative correlation function 

T{r, t) = (dt${x, t) ■ dt${x + r, t)) = {L/LfF{r/L) (5) 

where i^ is a new time-independent scaling function and 
L = dL/dt. Note that power-law growth, with or without 
additional logarithmic factors, implies that the prefactor 
[L/LY ~ 1/i^- If dynamic scahng holds both for T{r, t), 
as just defined, and for C{r,t), then the growth expo- 
nent can be determined through a self-consistent energy- 
scaling approach [ pTp^ . This restricted definition of dy- 
namical scaling, of both C{r,t) and T{r,tV picks up the 
scaling violations of the ID XY model |Q , and clearly 
separates the role of two-point from higher-point corre- 
lations [Q. We use this restricted definition here, and 
recommend it in the study of systems where dynamical 
scaling is questioned but Eqn. || seems to be satisfied. 



III. DYNAMICS 

We study purely dissipative quenches of 2D XY models 
from well above to below the Kosterlitz-Thouless tran- 
sition temperature Tkt- Because of the line of criti- 
cal points in the 2D XY model |12] the correlations in 
quenches to < T < Tkt have a modified scaling form 
p9| , p7[ . Essentially, critical equilibrium correlations have 
no characteristic length-scale and so the standard coarse- 
graining |l| to make temperature irrelevant to large-scale 
correlations is impossible. However, there is no indica- 
tion that temperature changes dynamical scaling, or its 
absence, in the 2D XY model. Accordingly, in this pa- 
per, we only investigate quenches to T = 0. The non- 
conserved coarse-grained dynamics [H are 



Fl 



d^x 



{V(i)f + Vo{4>^-if 



dt4> = -rSF/s^b, 

(<?(x,0)-^(x',0)) = A<5(x-x'), 



(6) 



where F is a kinetic coefficient that sets the time-scale, 
Vq is the potential strength that sets the 'hardness' of 
the vector spins, and A characterizes the initial disor- 
dered state. The orientation of the two-component order- 
parameter (/)(x) defines an angle 0(x) G [0, 27r], which is 
identical to the XY phase. The numerical implementa- 



tion of the dynamics is discussed below in Sec. IV A . 

In overview of the evolution: we start with a random 
high-temperature configuration and quench to T = 0. 
The order parameter locally equilibrates, but competi- 
tion between degenerate ground-states leads to topologi- 
cally stable vortices, with integer charges. The annihila- 
tion of oppositely charged vortices drives the subsequent 
dynamics, and characterizes one possible growing length 
scale — the vortex separation Ly. Of course, the order- 
parameter field around a moving vortex is not rigidly co- 
moving p7| ] , and so non-singular "spin- wave" distortions 
are generated by the dynamics even at T = 0. The dy- 
namics, emphasizing the vortices, can be visualized with 
a Schlieren pattern, see Fig. |l|, analogous to those used 
in the study of liquid-crystal films [E5[ . 



A. Scaling Correlations from Gaussian Closure 

Several approximation schemes eliminate high-order 
correlations in the evolution equation for two-point cor- 
relations |p], p5p^ , |35| , ^ . We use a Gaussian-closure ap- 
proximation, which gives quite good two-point correla- 
tions. We will use the results to normalize our corre- 
lations. This allows for a more sensitive test of scaling 
properties than has been possible before, and also high- 
lights weaknesses of this approach (see also |15 4^). 

For general 0{n) fields, we start with the Bray- 
Humayun-Toyoki (BPT) approach pa]. We introduce an 
auxiliary field rn parallel to the order parameter, fh — (/). 



The zeros of rh match the positions of the topological de- 
fect cores, while \rn\ is roughly the distance to the clos- 
est defect core. Assuming a Gaussian probability distri- 
bution for rh results in two-point correlations between 
(ri,ti) and (r2,t2): 



fl'y 

Cg{r,ti,t2) = — 



' 2' 2 



' 2' 2' 2 



-;7 



(7) 



where r = |r2 — ri|, B{x,y) is the beta function, 
and F{a, b; c; z) is the hyper-geometric function. The 
result is expressed in terms of the the normalized 
two-point, two-time correlation function of to: 7 = 
(to(1)to(2)}/[(to2(1)}(to2(2)}]1/2. 

The various approximation schemes schemes differ on 
the manner of determining 7. We use the the systematic 
approach introduced by Bray and Humayun |Q which 
produces 



"fir,ti,t2) 



Atih 



{h + hy 



d/4 



cxp(-rV[4(ti+i2)]), (8) 



where d is the spatial dimension. For equal-time corre- 
lations, we obtain the scaling form Cg{r,t) = fsprix), 
where x — r/L and L{t) = (4i)^/^. This highhghts a 
problem with all existing correlation-closure approaches 
as applied to 2D XY models, since while they recover 
a scaling form they miss the logarithmic factor in the 
growth-law p^ . The same scaling variable is used in the 
time-derivative correlation function 



Tgir,t) 



1 



16i2 



[j^x^Cy^ix) + j{x* - 4x2 ^ 2d)C^{x)] , 



(9) 



where C-y = dCg/dj and C^-y = d^Cg/dj"^. 



IV. SIMULATION 



A. Simulation Methods 



also [|13|,yj,gjj21|j2g]). The random initial conditions are 
chosen uniformly for each component from [—0.1,0.1]. 

We identify vortices with three methods that prove 
equally effective: by looking for the zeros in the vector 
field, by looking for plaquettes around which the phase 
rotates through ±27r, and by finding the peaks on the 
local energy density Ei = — J2\ 0(i) ' 4'ii); where the sum 
is over nearest neighbors of site i. Due to the periodic 
boundary conditions, the system has no net vorticity. 

In addition to tracking the number of vortices, we mea- 
sure several correlations of the "hardened" order param- 
eter, 0(j,t): 



C(r,i) = (0(j,i)-0(j+r,t) 



(11) 



The average (. . .} is over the independent sets of initial 
conditions, and includes a spherical average and an av- 
erage over lattice sites j. The structure factor is also 
calculated: 



Sik.t) = ^^{-k,t) ■ ^(k,t)j . (12) 

We also measure the time derivative correlation function. 



T{r,t)^(Stcf>{i)-St${i + r) 



(13) 



where Stfj) = </'(i + 1) ~ 0(0 i^ s- finite difference approxi- 
mation for the time derivative. 

To probe the distinction between vortex and non- 
vortex contributions to correlations, we measure a phase- 
gradient correlation function: 



= h{rlL)/L\ 



(14) 
(15) 



where the second line is the natural scaling ansatz for the 
correlations. Note that (V0) = 0. We then reconstruct 
the vortex contribution Dr {r, t) directly from the charges 
and locations of the vortices at a given time. From the 
vortex positions we build up the phase field V0(j) us- 
ing the periodic image of the minimal energy solution 
for each single vortex, V^9 — 0, due to Gr0nbech-Jensen 



We use a standard CDS update pS for soft spins, 
(/)(i, i), on a periodic lattice, where t is now a discrete 
integer time and i is the position: 

${h ^ + 1) = J I] [^(J. t) - ^ih t)] + E${i, t) tanh [|^(i, t) 
j 

(10) 

where cf) = <j)/ (f) is the unit vector. We use the stan- 
dard values D = 0.5 and E = 1.3. The dynamics are 
stable and have the same attractors as Eqn. a We do 
not observe pinning effects in quenches to T = (see 



dO / dx = — TT 



E 



sin(27ry)/[cosh(27r(x + n)) — cos(27r2/)]. 



d0/dy — n 






sin(27rx)/[cosh(27r(y + n)) — cos(27ra;)]. 



where x and y are the relative coordinates of the vortex 
in a system of size unity. The solutions for every vortex 
(with ±1 factors for vortices and anti-vortices, respec- 
tively) were added together for every point in the sys- 
tem to obtain the fully-periodic minimal- energy phase- 
field consistent with the vortex configuration. [Direct 
reconstruction of the order-parameter field (f> proved in- 
tractable due various counter-charge effects imposed by 



the periodic boundary conditions. In principle we could 
use our V0 reconstruction to recover the order-parameter 
field with additional line-integrations.] To obtain more 
accurate vortex positions, we first identify the lattice pla- 
quette by windings or energy peaks, then we use bilinear 
interpolation l^^l to more accurately locate the zero of 
the order parameter within the plaquette. The sign of 
the vortex is determined by the winding of the phase 
field around the plaquette. 



B. Simulation Results 

We simulate a size 512 x 512 system, averaging over 
40 independent samples. We check that there are no sig- 
nificant finite size effects in comparison to a 256 x 256 
size system, with 20 samples. The data for reconstructed 
correlations is currently restricted to the 256 x 256 size 
system. 

In Fig. H a), we plot C{r,t) with respect to distance 
scaled by L1/2 [C{Li/2,t) — 1/2]. The scaling is excel- 
lent, and the Gaussian-closure result {/bpt, solid line) 
is indistinguishable from the data. In Fig. b), however, 
the scaling collapse is not good with respect to the vor- 
tex separation [Ly, where p — l/L^] [Q. It is difficult to 
tell from this second plot alone whether scaling simply 
has a later onset time, or if scaling violations are indi- 
cated. This must be determined by a direct comparison 
of the length-scales L^ and ii/2, as well as by a study 
of the time-derivative correlations T{r, t), as discussed in 
Sec. 0. 

By normalizing the correlations with the Gaussian- 
closure result, Cg, we can sensitively probe scaling with 
the real-space correlations, see Fig. ^ While Cg is clearly 
too small at large scaled distance, correlations scale rel- 
atively well for i > 1000. 

The structure factor scales with respect to Lk = l/(fc}, 
its inverse first moment, see Fig. 0. Also shown (solid 
line) is the Fourier transform of the Gaussian-closure 
prediction, which slightly but systematically under and 
over-estimates the structure. By using a log-log plot we 
emphasize the S{k) ^ Pdefk~^ generalized Porod tail for 
k/{k) ^ 2, as per Eqn. ^. The good scaling of the Porod 
tail, which is determined by the vortex density, indicates 
that Lk ^ Ly asymptotically. 

We now directly test the assumption that all lengths 
asymptotically have the scaling growth-law of Eqn. g by 
plotting t/L'^ vs hit for L1/2, Lk, and Ly, in Fig. a. 
The scaling prediction is a linear plot, with non- universal 
slope and intercept given by the amplitude A and time- 
scale to ■ [Both of these can vary from one length-scale to 
another.] Linearity is observed for Int > 7.7 {t > 2200), 
in agreement with Fig. ^. We have fit them with straight 
lines with the same amplitude A but different to- The 
correlation length L1/2 has the strongest corrections to 
scaling, which is one cause of the bad scaling of C{r, t) 
when plotted vs. Ly in Fig. b). It is worth noting 



that the growing length scales can also be well fit using 
effective exponents of 0.42, 0.40, and 0.40 (±0.01), re- 
spectively, without logarithmic factors — see also |lq,p2| . 
However, if these effective exponents were asymptotically 
valid, and hence disagreed with the scaling prediction of 
Eqn. H we would not see scaling in the correlations [g7| . 

While the two-point correlations C{r,t) and S{k,t) 
support dynamical scaling, we must also investigate the 
time-derivative correlation function, T{r,t), as discussed 
in Sec. O. In Fig. P, we scale lengths with respect to 
Li/2, and remove the prefactor in Eqn. a by plotting 
r(r, t)/T{Li/2, t). While scaling only sets in for t > 2000, 
it is supported by the data. This correlation function has 
much more structure than the equal time correlations, 
such as a local maximum at x w 2.3 and a logarithmic 
divergence at small- a; due to fast vortex annihilations. As 
a result, it provides a more stringent test of the Gaussian- 
closure approximation. We find significant discrepancies, 
the first to be found in two-point correlation functions. 

Further confirmation of scaling in T{r, t) is found 
by exploring the time dependence of the amplitude 
T{Li/2,t) ~ t^^, see Fig. |^. The scaling form in Eqn.^ 
gives fj, — 2 (independent of logarithms), and we find 
/i = 2.0±0.1. This is consistent with scaling. In combina- 
tion with the scaling of C{r, t) [and S{k, t)], and the con- 
sistency of the growth laws of all measured length scales 
with the scaling result, we conclude that the quenched 
2D XY model dynamically scales. 

In the equilibrium 2D XY model, the singular (vortex) 
and non-singular (spin- wave) degrees of freedom have in- 
dependent contributions to the free energy |£l| . Could it 
be possible for such distinct 'singular' and 'non-singular' 
length-scales to exist in phase-ordering systems (see, e.g., 
[p2[)? If separation of vortices and spin- waves occurs, we 
expect spin-wave contributions to have a characteristic 
scale L ^ t^/'^, i.e. to have a faster decay with no loga- 
rithmic factor [E^. In which case, the direct correlations 
should either have scaling violations due to the differ- 
ent length scale or the spin-waves should be asymptot- 
ically irrelevant — leaving the direct and reconstructed 
correlations asymptotically equal at late times. As can 
be seen from the snapshots of |V6'| in Fig. g, the recon- 
struction maintains the vortex locations and is periodic. 
Indeed, the reconstruction provides the minimal energy 
configuration consistent with vortex positions — in other 
words any 'non-singular' contribution is absent. In Fig. ||, 
the correlation function for the direct and reconstructed 
fields are shown as a function of the scaled distance. We 
first notice that both correlations scale with respect to 
Ly — p^^i"^ but with different functional forms. Dr has 
a sharper knee at rp^" ~ 0.7, for example. This knee 
refiects the faster decay of V9 from the vortex core in 
the reconstructed configurations, as is apparent in Fig. ||. 
The significant differences between the bare and recon- 
structed correlations in the scaling limit indicate that 
both vortex and "spin-wave" contributions are relevant 
to the direct correlations, and that the separation seen 



in static properties does not hold in the dynamics. 

The Porod plot of the Fourier-transformed correla- 
tions, see Fig. ^, further highlights the differences (note 
the fc ^ 0+ intercepts). It is interesting that while 
(V^) ~ 0, the scaling curve has a non-conserved charac- 
ter. This is similar to correlations in globally-conserved 
systems. We also observe a k~^ Porod tail for kj p^l"^ > 2, 
which is expected from Fourier transforming the real- 
space scaling ansatz, Eqn. [l^, and setting the amplitude 
of the Porod tail proportional to the vortex density XjL? . 
The Porod tail has the same amplitude between the di- 
rect and reconstructed correlations, reflecting the singu- 
lar structure of the vortex cores. 



V. SUMMARY AND DISCUSSION 

We find no evidence for scaling violations in the 2D 
XY model. All lengths, L1/2, Lv, and Lk, have the 
same asymptotic form given by Eqn. 0, albeit with dif- 
ferent non-universal coefficients. Real-space correlations, 
structure, and time-derivative correlations all scale as ex- 
pected. Phase-gradient correlations, reconstructed from 
the vortex locations to have minimal energy and hence 
no spin- wave contributions, differed significantly from 
the direct correlations, indicating that both vortex and 
non-singular "spin-wave" contributions are asymptoti- 
cally relevant. We expect similar results to hold in closely 
related planar liquid-crystal systems. 

We have also shown how Gaussian-closure approx- 
imations can be useful to sensitively explore scaling. 
This has the added benefit of testing the approxima- 
tion schemes. In particular we find significant discrep- 
ancies with respect to the measured time-derivative cor- 
relations, T{r,t). More generally, we emphasize the role 
of sensitive null-like tests in checking apparent scaling 
violations. For example we plot the length-scales vs the 
expected growth law so that linear behavior is expected if 
scaling is obeyed. When scaling predictions are available, 
and in the face of transient corrections to scaling, this is 
preferable to the measurement of effective exponents. 

One can never absolutely rule out scaling violations, 
if only because simulations and experiments can never 
reach i — > 00. Each length in a system that dynami- 
cally scales will generically have different corrections to 
scaling. In quenches of the 2D XY model the leading 
correction is described well by toj the time-scale of the 
logarithmic factor. Since scaling violations seem to be 
rare in quenched systems, the assumption should be that 
systems dynamically scale without strong evidence to the 
contrary — including the inability to perform a scaling 
collapse with any length-scale for either C{r, t) or T{r, t). 
This provides a self-consistent confirmation of dynamical 
scaling, provided the lengths used for the collapse are 
consistent with the same asymptotic growth-law. For 
some systems, including this one, the scaling growth- 
law can be independently determined. This is invaluable 



when long-lived (logarithmically decaying) corrections to 
scaling are expected. The scaling of some other lengths 
in the problem can sometimes also be required for con- 
sistency. In this case the defect-separation scale Ly is 
needed to set the Porod amplitude, and hence must be 
consistent with the lengths L^ and ii/2 extracted from 
S{k,t) and C{r,t), respectively. 
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FIG. 1. Schlieren patterns in various times after a quench 
of the 2D XY model in a size 128 x 128 system. The inten- 
sity is sin^(2^), where is the local XY phase. Each vortex 
emanates 8 brushes, alternating white and black. 




FIG. 2. a) Real-space correlations vs a; = r/Li/2, where 
C{Li/2) = 1/2. The continuous curve represents the theoret- 
ical prediction, fsPTix). b) Attempted scaling with respect 
to the vortex density. 
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FIG. 3. The difFerence between measured correlations, 
C{r,t), and the Gaussian-closure prediction, Cg, normalized 
by Cg and plotted against scaled distance. 
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FIG. 5. We plot t/L^ vs Ini for three lengths: I/1/2, ifc, 
and Lv The observed linear dependence at late times indi- 
cates that the dynamical scaling growth law, Eqn. 0, holds. 
As shown by the parallel straight-lines, the offset (given by 
to) is non-universal. 
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FIG. 4. The structure factor in a log-log Porod plot. The 
first moment, (k), is used to rescale momenta. The continu- 
ous line is the Gaussian-closure prediction. Symbols are the 
same as the previous figure. 
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FIG. 6. Scaling plot of T{r,t)/T{Li/2,t) vs r/Lj/a. The 
continuous curve is the theoretical prediction of the Gaussian 
closure scheme — significant deviations are apparent. 




FIG. 7. Log-log plot of the time dependent prefactor of 
the time-derivative correlation function. The best fit over the 
range shown yields a decay 1/t^ with an exponent ^ = 1.96. 
Varying the fit range yields /i = 2.0 ± 0.1. 
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FIG. 9. Direct and reconstructed V6 correlations, D{r,t), 
and Dr{r,t), respectively, normalized by vortex density for a 
scaling plot vs scaled distance. Scaling is observed after the 
earliest time. 
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FIG. 8. Snapshots of \V0\ for a 256 x 256 size system at 
two times after the quench. The left column shows the direct 
|V6l| , with contour levels at 0.1, 0.2, 0.4, 0.8, and 1.6. [Note 
that the lattice spacing defines a unit of length, so the largest 
gradient magnitude is tt.] The right column shows |V6'| pe- 
riodically reconstructed using only the vortex positions, with 
the same contour levels. Significant differences between the 
direct and reconstructed VS field can be seen away from the 
immediate vicinity of the vortices. 
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FIG. 10. Porod plots of Fourier-transformed direct and 
reconstructed V6 correlations, D{k,t) and Dr{k,t), respec- 
10 lively. Symbols are the same as the previous figure. 



